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1 Introduction 

In classical mechanics Hamiltonian system on a 2n-dimensional phase space M is called com- 
pletely integrable in Liouville's sense if it possesses n functionally independent integrals of motion 
Hi, . . . , Hn in involution: 

J IT 

-^ = {H,Hi} = 0, {Hi,Hj} = 0, i,j = l,...,n, 

where H = Hi is the Hamilton function and {•, •} is the Poisson bracket on M. 

Superintegrable system is a system that is integrable in the Liouville sense and that possesses 
more functionally independent integrals of motion than degrees of freedom. The construction of 
superintegrable Stackel systems using angle variables Uk has been proposed in [12, 13, 14, 15]. 

In generic case the action variables cj^ are multi- valued functions on the whole phase space M. 
In fact, we can extract polynomial integrals of motion from angle variables only when we can 
apply addition theorems to the corresponding Abelian integrals. As there are only few addition 
theorems for the Abel equations [1, 3] we can easily classify the corresponding superintegrable 
systems, see [12, 13, 14, 15]. 

The goal of this brief note is to present some trivial examples of applying this generic theory 
associated with elliptic and parabolic coordinate systems on the plane. Superintegrable systems 
separable in spherical coordinates can be found in [2, 6]. The corresponding addition integrals 
of motion are related with an addition theorem for the logarithmic angle variables. Of course, 
there is a trivial generalization of the proposed method for all the orthogonal coordinate systems 
in (ellipsoidal, paraboloidal, cylindrical, prolate and oblate spheroidal coordinates etc). 

The non-Stackel superintegrable systems in classical and quantum mechanics have been con- 
sidered in [7, 9] . In contrast with the Stackel case we do not have a generic theory for constructing 
such superintegrable systems. 



*This paper is a contribution to the Special Issue "Superintegrability, Exact Solvability, and Special Functions" . 
The full collection is available at http://www.emis.de/journals/SIGMA/SESSF2012.html 
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2 The Stackel systems 

The system associated with the name of Stackel [11] is a holonomic system on the phase 
space M^", with the canonical variables q = (qi, ■ ■ ■ , qn) and p = {pi, . . . ,pn)- 

n 

J7 = ^ dpj A dqj , {pj , qk} = 5jk. 
j=l 

The nondegenerate n x n Stackel matrix S, whose j column depends only on coordinate qj, 
defines n functionally independent integrals of motion 

n 

Hk = Y.iS-')^,ip] + V,{q,)), 
i=i 

the separated relations 

n 

p'j=Yl HkSkjiqj) - Vjiqj), 
k=l 

and the action variables Wk 

1 ^ SijW 



n 

i - 



n 



E/ I =d^- (2-1) 



lj:HkSkj{X)-Vj{X) 

k=l 



SO that 

{Hj, Hk} = {u}i,ujk} = 0, {Hj,LOk} = Sij. 

In generic case the action variables (2.1) are sums of the multi-valued Abelian integrals. 
However, if we are able to apply the known addition theorems for these Abelian integrals then 
we can get additional integrals of motion [12, 13, 14]. 

Let us discuss addition theorems for the logarithmic (exponential) and elliptic functions [1, 3]. 
In the first case polynomials 

n 

P, = Y,HkSkj{X)-Vj{X) 

k=l 

are the second-order polynomials and 



1 

2^ I U \ nk^ 



1 ^ dX ( , kjq^+h, 



kn^ + bjX + cj ,=1 



■1 



So, we can easily get polynomial or rational function 



at the special choice of kj and z. In the second case m — 1 angle variables w^, . . . , Wfc+m-i with 
1 < m < n have to satisfy to the Abel equations 

dxi dx2 dxm , 

H ]== H 1 , , = diOk, 
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+ , . . + ••• + 



+1) 



m-2 
^1 



dxi x^ ^dx; 



,m-2 



+ 



+ ••• + 



dx^ 



fc+m— 1) 



with a common polynomial of fixed degree 2m 

fix) = ^2mx2" + ^2n^-lx2'"-l + • • • + ^iX + ^q- 

If this case there are some additional Richelot integrals of motion [10] 



(2.3) 



Ck 



F'(xi) afc-xi 



+ 



+ 



Clfc x^ 



\//(xi) Vfi^rn) 



F'(xi) 



2m 



Here are values of x at the branch points of the corresponding hyperelliptic curve and F{x) = 
(x - xi)(x - X2) • • • (x - Xm) [10]. 

At m = 2 we have famous Euler algebraic integral [3]. If = and A2m-i 7^ in (2.3) 
there is another additional Richelot integrals of motion [10]. The Weierstrass generating function 
of such integrals for any values of the coefficients and other constructions of additional 
integrals of the Abel equations are discussed in [1]. Some example of the Euler and Richelot 
superintegrable systems may be found in [4, 12, 15]. 

Below we show how these addition theorems could help us to classify superintegrable systems. 



3 Elliptic coordinate system 

Let us consider elliptic coordinates on the plane gi,2 defined by 



1 



X 



X 



^ (A-gi)(A-g2) 
X + K A2 - k2 



K G 



The corresponding momenta reads as 

2p^x 2pyy 
Pi = \ P2 



2pxX ^ 2pyy 



qi - K qi + k' ' " q2 - K q2 + n 
The Stackel matrix and the separated relations 



S 



ql 



1 



g2 \ 

9 

92 - 
1 



Pj + Vl 
PI + V2 



qiH 



ql 



q2Hi 



give rise to the following Hamiltonians in the involution 
Hi 



{ql-^''){pl + Vi) . {ql-^''){pl + V2) 

I 



q2 - qi 

,2 ,^2\f„2 



q2 - qi 

2 ^2\f„2 



Ho 



q2{qi - K2)(pf + Vi) qi{ql - K^){pl + V2) 



qi - q2 



qi - q2 



H2 



q\ 



Ho 



ql 
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The Hamiltonian Hi commutes with the second angle variable W2, which is equal to 
1 d\ 

W2 = 



1 dX 

Polynomials 

Pi 2 = (A^ - K^) {XHi + H2- yi,2A2 + ^1,2^^) (3.2) 

standing under square root in these integrals are at least third-order polynomials on A. So, in 
this case we can not apply addition theorem for the logarithms. 

It is easy to see that we can apply addition theorem for the elliptic functions at 

Vi = Vi = a. 

Namely, if we put A = x and A = y in the first and second integrals (3.1), we could apply the 
Euler addition theorem 

dx dy ds , , 

^ + ^ = -= 3.3 



to angle variable 012 ■ Here X is an arbitrary quartic 

X = ax^ + 46x3 + 6cx2 + idx + e (3.4) 

and Y, S are the same functions of another variables y and s. In this case, symmetrical bi- 
quadratic form of X and y 

F(x, y) = ax^y^ + 26xy(x + y) + c(x^ + 4xy + y^) + 2d(x + y) + e = 

defines the conic section on the plane (x, y). According to [1, 3, 10], there is a famous Euler 
integral 

F{x,y)-VxVY i/VX-VyV a(x + y)^ 

C = o/„ = 7 „ . 6(x + y)-c. (3.5) 




For the quartic (3.2) associated with the angle variable (3.1) this Euler integral looks like 

rr (.Pi - P2)(.qi - K^){q2 - l^"^) f ( .2 ( \2 2 ^ ( ^ 

Hz = 7 To a{qi - q2) - [pi - P2 ^ + [Piqi - ^2^2 

[qi - q2y V 

Here qi^2 and pi^2 are elliptic coordinates and momenta. 

It is a third-order polynomial in momenta which commutes with the Hamiltonian 

{/7i,i72} = 0, {H2,H;} = H^^Q. 

The algebra of the polynomial integrals of motion Hi, H2, H^ can be closed only after some 
other polynomial generators are added. 

Thus, we easily find the additional integrals of motion for the Hamilton function of the 
oscillator 

iHi = pI + pI + a{x^ + y^) 

using the separation of variables in elliptic coordinate system and the corresponding angle va- 
riables. Another result is that there is only one superintegrable system separable in elliptic 
coordinates and associated with the known addition theorems for Abelian integrals. 



Superintegrable Stackel Systems on the Plane: Elliptic and Parabolic Coordinates 



5 



4 Parabolic coordinate system 

Let us consider parabolic coordinates on the plane qi^2 defined by 

X = qm, y = — ^ — 
and the corresponding momenta 

_ Pi 92 + P2qi _ qipi - q2V2 

2,2' Py 2,2' 

qt + q2 qi+ qi 

The Stackel matrix and the separated relations 

g^fql ql\ pl + VM-qlH,-H2 = 0, 

VI -V' pi + 1/2(92) -92^1 +^2 = 0, 



give rise to the Hamiltonians 

r / _ \ __2_2 __2_2, _2ta/_ n '^''^ ^ J2j 



Pi+ptViiqi) + V2iq2) rr _ Pi92-P29i+ 92^1(91) -91^2(92 



2 I 2 ' 2 1 2 

91 + 92 91 + 92 
The Hamiltonian Hi commutes with the second angle variable W2, which is equal to 
_ 1 dX 1 dX 

^^^"27 ^x2Hi + H2 ^Wa) y/xmi-H2^V2W' 

In contrast with the elliptic coordinates, polynomials 

Pi,2 = A^Fi + F2 -yi,2(A) 

standing under square root in these integrals are at least second-order polynomials on A. So, we 
can apply both known addition theorems to these Abelian integrals. 
In fact, these integrals are expressed via logarithmic functions iff: 

Case 1 : Vi = biqi + ci, V2 = 6292 + 02, 
Case 2 : Vi = aigf^ + 6i, ^2 = '2292'^ + ^i- 
The addition theorem for the elliptic function is applicable iff: 

Case 3 : Vi = aiqi + biq^ + cigf , V2 = ai92 ~ ^i92 + ^192" • 

The corresponding Hamilton functions are deformations of the Kepler-Coulomb and oscillator 
Hamiltonians: 

Case 1 : Hi=pl+pl^ , ^ ( 6iv/x + + + 62 a/x/x^ + ^ + ci + C2 | 

2Jx^ + V / 



2r\ 2 

Case 2 : Hi=pl+ pi + — ^ \ { + +61 + 62 



pI+pI + -^ (6ix/2cos^ + 62V2sin^ +ci + C2 



2\/ x"^ + y"^ I X + -\/x2 + x — -s/ x'^ -\- 



2 



^2,2, 0^ <^ 61 + 62 

Px +Py+ (cos (/?+!) 2r2 (cos ^ - 1) 2r ' 

Case 3 : Hi=p1+pl + a (4x^ + y^) + 2/3x + ^. 

Here r = y^x^ + y2 ^nd 99 = arctanx/y are polar coordinates on the plane. According to [8] 
these systems remain superintegrable in the quantum case. 
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4.1 Case 1 

In the first case the second angle variable equals 



W2 



<?i 



dX 



1 



VX'^Hi -biX + H2-ci 2 J ^X^Hi - 62A - H- 



12 



dX 



C2 



1- (pi - '^^^^) In (p2 - '^^^) 



The application of the addition theorem (2.2) to 002 gives rise to the following rational integral 
of motion 

Z = e2Vlhco2 ^ 2qiHi + 2pi^i - bi 
2q2Hi+2p2y/H[-h2 

In order to calculate polynomial integral of motion let us consider a series expansion of the 
function 



/ 



1 . y^oy-i-, am,H2-mci-bi)-f3i4H^H2 + 4H,C2 + bl) , 



'Hi 



^Hi{AHiC2 + 4H1H2 + bj){b'i + AHiCi - AH1H2) 



by momenta |?i,2- Here we substitute the variables qi^2 from the separated relations (4.1) into 
the rational integral Z and a, (3 are undefined polynomials in Hi 2- 

Equating first coefficient of this expansion to zero one gets the following expressions for these 
polynomials 

a = AH1H2 + 4HiC2 + bl /? = 4H1H2 - AHiCi - bj. 
At this values of a and f3 the function / becomes a third-order polynomial in momenta 



^3 



1 



1 



ae 



+ (3e 



8(pig2 -P2gi)(Pl +P2 +Cl +C2) , 4(2^151^2-^2(91-92))^! 



ql + ql 



+ 



H2p2qiq2 + Piiqf - q2))h 



2 1 2 ' 

qt + Qi 

such that 

{Hi,Hs} = 0. 

In order to close the algebra of the polynomial integrals of motion Hi, H2, H3 we have to add 
one more polynomial generator 

H4 = {H2, Hs} = 2ae2v^'^2 - 2^e-^^'^\ 

by analogy with exp(a;), sin(a;) and cos(a;) functions. 

Remark 1. One of the referees proposed another construction of the polynomial integrals of 
motion from the angle variable lj2. Namely, from the separation equations we can deduce that 



Ai = i2qiHi + 2pi - bi){2qiHi - 2pi - 61) 
A2 = (2g2^i + 2p2^i - b2)i2q2Hi - 2p2^i - 62) 



-4HiH2+iHiCi+bl 
4H1H2 + 4HiC2 + bl 
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We can therefore write 



and 



^-1 = AiZ-i = {{2qiHi - h){2q2H2 - 62) - 4Hmp2) 

- 2^/lh{{2q2Hi - b2)pi - {2qiHi - b{)p2) 



^2 = Aa^ = [{2qiHi - 6i)(2(72^2 - ^2) - ^Hipip2) 

+ 2^i{{2q2Hi - b2)pi - {2qiHi - bi)p2). 

Consequently 

^3 = -^{^1 - ^2) = -4((2g2i^i - b2)pi - {2qiHi - 61)^2) 

is a third-order constant of motion. This method of explanation may be clearer than the method 
of expansion of the rational in momenta function / = aZ + j3Z~^ with indefinite coefficients a 
and /3. 

Remark 2. Let us remind, that two-dimensional open Toda lattice defined by the following 
polynomial integrals of motion 

^1 =P?+P2 + e^^"''^ H2=Pl+P2, 

has the non-rational in momenta additional integral of motion 



Z = ^^-^^ + ^ expf^/j^^V J = 2H2-Hl 

P1-P2-VJ V P1+P2J 



Pi -P2 

which can be also obtained from the second angle variable U2 [14]. However, it is easy to prove, 
that we can not apply working above constructions of polynomial integral of motion H3 in this 
case. 

4.2 Case 2 

In the second case the angle variable is equal to 

1 XdX 1 XdX 

UJ2 



27 y'X^Hi + {H2 - 6i)A2 - ai 2j ^ X^Ri - {H2 + b2)X^ - a2 

Changing variables // = one gets second-order polynomials 

Pj = iJ?Hi ± {H2 =F bj)ii - Qj 

under the square root and desired sum of the logarithms 

In {2qfHi + 2qipi^i + H2 - 61)) In {2qlHi + 2q2P2VH~i - H2 - 62)) 
L02 = 



The rational integral of motion (2.2) is equal to 
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As above we consider the expansion of the function 

/ = -^(aZ + /3Z-i) 
V-til 

by momenta pi,2- Equating first coefficient of this expansion to zero one gets polynomials a, (3 

a = AHia2 + + 262^2 + H^, (3 = -AHiai -b\ + 2hH2 - H^. 
At this values of a and /? the function / becomes a third-order polynomial in momenta 

H3 = {aZ + PZ'^) = -L (ae^v^^2 + jie'^^'^A 

4^{qiPi - q2P2)iq2Pi - qiP2f , 4^qiq2iq2Pi - qiP2){bi + b2) 



+ 



qf + q2 qi + ql 

^ 4aig2(g2giPi - + ^2)^2) _ 4a2gi(g2giP2 - (^gj + 

ql{ql + ql) ql{ql + ql) 

such that 

{/7i,i73} = 0. 

In order to close the algebra of the polynomial integrals of motion Hi, H2, H3 we have to add 
one more polynomial generator 

Hi = {H2, H3} = 4ae^^'^2 _ 4/3e-4v^^2_ 
4.3 Case 3 

In the third case the angle variable is equal to 

XdX XdX 



bJ2 



XdX r 

J V-aiX^ - 61 A6 + HiX^ + H2X'^ -ci~ J V-aiX^ + 61 + HiX'^ - H2X'^ - ci ' 



Changing variables A = \/x and A = i^/y at the first and second integral one gets the Euler 
addition theorem (3.3). In fact, this example has been considered in Euler's book [3] too. 
Identifying quartic 

P = -oi/i^ - 61^^ + Hin^ + H2ti - ci 

with X (3.4) we can easily calculate the Euler integral of motion (3.5) in parabolic coordinates 

^ ^ {qipi - q2P2){qiP2 + 92^1)^ aigig2(2gfy>2 + q2qlpi - qiqlp2 - 2g|pi) 
' ' {ql + ql? ql + ql 

biqiq2{qiP2 + q2Pi) ci{qipi - 92^2) 

~r 0.0 ' 



qf + qi qlql{ql + qD 

The algebra of the integrals of motion Hi, H2, H^ is more complicated then the algebra asso- 
ciated with the addition theorem for logarithms. In fact, in order to close this algebra we have 
to introduce the counterparts of the Jacobi elliptic functions sn(a;), cn(a;) and dn(a;) instead of 
the trigonometric functions sin(a;) and cos(a;), which we used for the superintegrable systems 
associated with the addition theorem for logarithms. 
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5 Conclusion 

It is known that orthogonal coordinate systems on Riemaniann manifolds can be viewed as an 
orthogonal sum of certain basic coordinate systems and these basic systems can be obtained from 
the elliptic coordinate system [5] using a degeneration procedure. This degeneration decreases 
the degree of polynomials standing under square roots into the angle variables (2.1). Thus, we 
have only one superintegrable systems separable in elliptic coordinates, whereas for degenerations 
we have a lot of different superintegrable systems. As usual, the addition theorem for logarithms 
allows us to get additional integrals of higher order in momenta [12, 13]. 
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